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We propose an efficient broadband frequency generation technique for two collinear optical para- 
metric processes ljs — Ui + LJ2 and 0)4 = u\ — u)2- It exploits chirped quasi-phase- matched gratings, 
which in the undepleted pump approximation regime perform population transfer that is analogous 
to adiabatic population transfer in a three-state "vee" quantum system. The energy of the input 
fields is transferred adiabatically either into ljs or C1J4 field, depending on which of the two phase 
matchings occurs first by the local modulation period in the crystal. One can switch the output 
between uj$ and 0J4 by inverting the direction of the local modulation sweep, which corresponds to 
a rotation of the crystal by angle ir. 

PACS numbers: 42.65.-k, 42.65.Ky, 42.79.Nv, 42.25.Kb 



I. INTRODUCTION 

Sum frequency generation (SFG) and difference fre- 
quency generation (DFG) occur, when two input beams 
generate another beam with the sum (SFG) or the differ- 
ence (DFG) of the optical frequencies of the input beams 
[lHH . In order to be efficient these processes traditionally 
require phase matching , which is usually difficult to 
achieve simultaneously for SFG and DFG [lpy]. 

Recently Suchowski et al. [IHI] used an aperiodically 
poled quasi-phase-matching (QPM) crystal to achieve 
both high efficiency and large bandwidth in SFG and 
DFG in the regime of undepleted pump approximation. 
Their approach was based on ideas from rapid adiabatic 
passage in two-state quantum systems [7H9(. 

In this paper, this method is further extended to re- 
alize a potentially highly efficient broadband SFG and 
DFG with a single crystal. This is achieved by treating 
the two simultaneous collinear second-order parametric 
processes u> 3 = Wi + uj 2 and lq± = u\ — uj 2 in analogy to 
coherent population transfer in three-state "vee" quan- 
tum systems. The analogy to level crossings in atomic 
systems ITHSl ensues from the linearly chirped QPM grat- 
ings [HHL2I ] that we use to this end. 



II. BACKGROUND 



id z E 1 = Oi£;*£;4e" iA2Z , (2a) 
id z E 2 = n 2 ElE 4 e- lA2Z , (2b) 
id z E 4 = Q,iEiE 2 e iA2Z , (2c) 

where z is the position along the propagation axis, c is 
the speed of light in vacuum, and Ej, u>j and rij are the 
electric field, the frequency and the refractive index of the 
j-th laser beam, respectively. The coupling coefficients 
are 

% = x {2} 0J 3 /Acn 3 (j = 1,2,3,4), (3) 

while the phase mismatches for SFG and DFG processes 
are 

Ai = U\ri\/c + u 2 n 2 /c - uj 3 n 3 /c + 2n/A, (4a) 
A2 = n\U)\jc + uj 2 n 2 /c — W4714/C + 2ir / A. (4b) 

We combine Eqs. (JTJ) with Eqs. ^ to write a system 
of differential equations that describes the two simulta- 



neously running processes 

id z E x = £li (ElE 2 e- tAlZ + ElE A e~ lA2Z ) , (5a) 

id z E 2 = n 2 {E{E z e^ z + ElE 4 e- tA2Z ) , (5b) 

id z E 3 = n^E^^, (5c) 

id z E 4 = 4 £i£ 2 e ,A22 . (5d) 



The SFG and DFG processes for a QPM crystal with 
susceptibility an d local modulation period A(z) are 
described by two sets of nonlinear differential equations 

m 

id z E x = Q, x E* 2 E z e- iAlZ , (la) 

id z E 2 - ft 2 £*£ 3 e- lAlZ , (lb) 

id z E 3 = Q 3 E 1 E 2 e iAlZ . (lc) 



III. ADIABATIC EVOLUTION OF SFG AND 
DFG IN UNDEPLETED PUMP 
APPROXIMATION 

The coupled nonlinear equations ([5]) can be linearized 
if we assume that the incident pump field E\ is much 
stronger than the other fields. In this case its amplitude 
is nearly constant (undepleted) during evolution and as 
a result Eqs. ([5]) arc reduced to a system of three linear 
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equations, 



id z A(z) = H(z)A(z), 
Ai Qp 

n P o 
o n s 



o 

A 2 



(6a) 
(6b) 



with 



n P 


— Ei^/n 2 n 3 , 


(7a) 




= Ei \j r^2^4, 


(7b) 


A(z) 


= [A 3 (z),A 2 (z),A i (z)] T , 


(7c) 




= E 2 /^n~ 2 , 


(7d) 


A, 


= Eze-^/^h, 


(7e) 


A 4 


= E 4 e- iA2Z /-s/Sh. 


(7f) 



Upon the substitution z — > t, Eq. ([5]) becomes identical 
to the time-dependent Schrodingcr equation for a three- 
state quantum system in the rotating-wave approxima- 
tion; the vector A(z) and the driving matrix H corre- 
spond to the quantum state vector and the Hamilto- 
nian, respectively. The diagonal terms of the matrix 
H, Ai, and A 2 , correspond to the detunings while 
the off-diagonal terms, flp and Qs, correspond to pump 
and Stokes Rabi frequencies. We note that the quan- 
tity |A(z)| 2 = \A 3 (z)\ 2 + \A 2 (z)\ 2 + \A i (z)\ 2 is conserved, 
analogously to the total population in a coherently driven 
quantum system. If the energy is initially in the input 
electric field with frequency uj 2 



A= [0,A 2 ,0] 



(8) 



then the three linear equations © will form a "vee" pat- 
tern analogously to the "vee" configuration of a three- 
state quantum system (see Fig. [I). We assume that 
phase mismatches either increase (sign +) or decrease 
(sign — ) linearly along z 



Ai = 5i ± a 2 z, 
A 2 = S 2 ± a 2 z, 



(9a) 
(9b) 



which can be achieved, for example, by varying the lo- 
cal modulation period A(z). For the sake of generality, 
we take hereafter a as the unit of coupling and 1/a as 
the unit of length. Therefore the three eigenvalues of H 
cross at two different distances z m (m = 1,2), thereby 
creating a crossing pattern in analogy with two paral- 
lel energies crossed by a third, tilted energy in quantum 
physics [l3j . This crossing pattern can be easily exam- 
ined by the famous Landau-Zener-Stiickelberg-Majorana 
(LZSM) model EM3, which is the most popular tool 
for estimating the transition probability between two 
states of crossing energies. This model assumes a con- 
stant interaction of infinite duration and linearly evolv- 
ing energies. Owing to some mathematical subtleties, the 
LZSM model often provides more accurate results than 
anticipated. Its popularity is further promoted by the 




FIG. 1: (Color online) The "vee" linkage pattern for the lin- 
earized system ((6J, treated in analogy to the "vee" configura- 
tion of a three-state quantum system. 



extreme simplicity of the transition probability expres- 
sions it derives. The LZSM model has been extended to 
three and more levels by a number of authors. In the 
Demkov-Osherov (DO) model [3 (H, a single tilted en- 
ergy crosses a set of N parallel energies. Our case, cf. 
Eqs. ([HI) and Eqs. ©, matches the DO model with two 
parallel energies crossed by single tilted energy, as shown 
in the top frames of Fig. [5J 

The proposed SFG and DFG are illustrated in Fig. H 
The top frames plot the eigenvalues of H (ei, e 2 and £3) 
vs z. Initially only the uj 2 field is present (see Eq.©). 
If the evolution is adiabatic then there are two possible 
paths that the system can follow (left and right frames). 
If the phase match for the uj 3 generation process occurs 
first (left frames of Fig. [2]), then the energy is passed into 
the W3 field. If instead first we observe the phase match 
for the L04 generation process (right frames of Fig. [lj, 
then efficient energy transfer to the W4 field takes place. 

The bottom frames of Fig. [2] show the evolution of 
the normalized light intensities for the three possible fre- 
quencies w 2 , w 3 and W4. The left and the right scenarios 
extend the single-step adiabatic passage scenario for ei- 
ther SFG or DFG 0Q. 

In the beginning and at the end each eigenfrequency 
£i(z), i = 1,2, 3, coincides with one of the intrinsic diag- 
onal terms of H, while in between it is a superposition 
of these terms. In the adiabatic limit, the system follows 
the eigenstate of H(z), which asymptotically coincides 
with the initial state of the system. Correspondingly, 
the frequency of the system at any instant of z is the 
frequency of this state, i.e. for the initial condition © 
E\(z). However, the composition of £1(2) is different for 
DFG and SFG because the order of the crossings differ. 
Reordering the crossings can be easily done by z reversal. 
Hence one can achieve either SFG or DFG with a single 
chirped QPM crystal just by rotating it on an angle tt. 

Next we turn our attention to the conditions needed 
for adiabatic evolution in the two distinct cases of SFG 
and DFG. By applying the LZSM model [TM^] 



p = 1 - exp(-27rft 2 /" 2 ), 



(10) 



we find that to obtain transition probability larger than 
1 — e we must satisfy the following conditions at each 
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FIG. 2: (Color online) Schematic evolution of frequency gen- 
eration obtained by numerical integration of Eqs. ([6]) for 
Si = —10a, 82 = 5a and Q p — Q, s = 2a. Left frames: SFG 
with increasing phase mismatches. Right frames: DFG with 
decreasing phase mismatches. Top frames: Diagonal elements 
(solid lines) and eigenvalues (dashed lines) of the driving ma- 
trix H of Eq. Bottom frames: normalized light intensi- 
ties. The sign difference in Eq. © between the left and the 
right frames corresponds to z reversal (rotation of the crystal 
by angle ir). 



readily verify that the conditions arc satisfied for the pa- 
rameters used in Fig. [2] 

As was shown earlier [H-Q m analogy with atomic 
physics, adiabatic implementation of both SFG and DFG 
leads to robustness of the adiabatic approach against 
variations of the parameters such as propagation dis- 
tance, couplings and initial (final) phase mismatches. 
Our approach shares the same robustness as previously 
proposed adiabatic schemes, which includes stability to 
variation of crystal temperature, wavelengths of the in- 
put electric fields, crystal length and angle of incidence. 

IV. CONCLUSION 

We have used the analogy between the time-dependent 
Schrodinger equation and the SFG /DFG equations in the 
undcpleted pump approximation regime to propose an 
efficient broadband SFG/DFG technique realized with a 
single crystal. A local modulation period sweep along the 
light propagation creates crossings in the phase matching 
between different parametric processes, which in combi- 
nation with adiabatic evolution conditions allow for both 
efficient and robust SFG and DFG for the input frequen- 
cies. Chirped adiabatic QPM gratings offer robustness 
against variations of the parameters of both the crystal 
and the electric fields, which include the crystal temper- 
ature, the wavelengths of the input electric fields, the 
crystal length and the angle of incidence. 

The present work can be viewed as a generalization of 
the idea of Suchowski et al. 0, however applied simul- 
taneously for SFG and DFG in a single crystal. 



crossing 



> 



a 



ln(l/e) 



2- 



(11) 



where Q, = Q, p for SFG and fl = fls for DFG. One can 
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